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Abstract 

We study the collective excitation spectrum of a d = 3 site-disordered 
Anderson-Hubbard model at half-filling, via a random phase approximation 
(RPA) about broken-symmetry, inhomogeneous unrestricted Hartree-Fock 
(UHF) ground states. We focus in particular on the density and charac- 
ter of low-frequency collective excitations in the transverse spin channel. In 
the absence of disorder, these are found to be spin wave-like for all but very 
weak interaction strengths, extending down to zero frequency and separated 
from a Stoner-like band, to which there is a gap. With disorder present, a 
prominent spin wave-like band is found to persist over a wide region of the 
disorder-interaction phase plane in which the mean-field ground state is a dis- 
ordered antiferromagnet, despite the closure of the UHF single-particle gap. 
Site-resolution of the RPA excitations leads to a microscopic rationalization 
of the evolution of the spectrum with disorder and interaction strength, and 
enables the observed localization properties to be interpreted in terms of the 
fraction of strong local moments and their site-differential distribution. 
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I. INTRODUCTION 



Central to the understanding of any system with long-ranged order is a description of 
its collective excitations. If the Hamiltonian has a continuous symmetry which is broken in 
the ground state, their spectrum is gapless (as implied by the existence of Goldstone modes) 
and they dominate low-temperature thermodynamic properties. 

One of the most familiar examples of such a system is the Heisenberg model, whose 
collective excitations are pure spin waves. With S = |, antiferromagnetic (AF) nearest- 
neighbour interactions and on a bipartite lattice, the ground state of the Hamiltonian has 
long-ranged magnetic order in three dimensions M. While no exact solution exists for d > 1, 
a successful and longstanding approach to the ground state and collective excitations of the 
model is given by linear spin-wave theory (LSW)i, in which particle-hole pairs created out 
of the Neel mean-field ground state are treated as bosons. The theory becomes exact in the 
limit d —>■ oo, and 1/d corrections can be systematically incorporated!^ but even in d — 2, 
LSW estimates for the sublattice magnetization, spin wave velocity and ground state energy 
are in good agreement with quantum Monte Carlo calculations (QMC) (see e.g. Ref. |3] for 
a comprehensive review). 

It is well known that the AF Heisenberg model represents the strong coupling limit of the 
half-filled Hubbard model!. At finite interaction strengths U, the mean-field ground state 
remains the Neel state (albeit with self-consistently determined local moment magnitudes) 
and gives a qualitatively sound description of the true ground state in d > 2. The collective 
excitation spectrum is however considerably more complex than that of the Heisenberg 
model, since its low-frequency spin wave excitations are mixed with incoherent, Stoner-like 
processes, leading to 0(N 2 ) particle- hole excitations in contrast to the N spin waves which 
alone survive in the strong-coupling limit. 

The generalization of LSW to finite interaction strengths is given by the random phase 
approximation (RPA) about the fully unrestricted Hartree-Fock (UHF) saddle point. Al- 
though naively thought of as a weak-coupling theory, it reduces precisely to LSW in the 
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strong coupling limit, an observation which has led a number of authors to conclude it to 
be a sound way of probing the collective excitations of the Hubbard model over a wide 
range of interaction strengths: for example, such an approach has been successfully used to 
investigate the transverse spin excitation spectrum of the Hubbard model both in d = li 
and d = 2 (see e.g. Refs [7], []). 

In the present work we consider initially the collective excitation spectrum of the Hubbard 
model on a simple cubic lattice. This is shown to possess two distinct bands for all interaction 
strengths U: a low-energy spin wave-like band corresponding primarily to orientational 
fluctuations of the local moments, extending down to zero frequency and containing two 
Goldstone modes; and an upper band of Stoner-like excitations, to which there is a gap due 
to Fermi-surface nesting in the corresponding single-particle spectrum0. The latter involve 
a significant degree of charge transfer character (site double-occupancy) and are projected 
out in the large U limit, leaving solely the low-frequency excitations, which become the pure 
spin waves of the Heisenberg model. It is found that, despite mixing between the two bands 
at finite U, the low-frequency excitations retain a strong degree of spin wave character down 
to weak interaction strengths. 

A problem about which very little is known, even at a qualitative level, is the effect 
of disorder on the collective excitation spectrum. In a disordered system, the collective 
excitations are governed by the interplay of a number of complex processes, including strong 
correlations of particle-hole pairs, Anderson localization of the underlying single-particle 
states, and the non-trivial statistics associated with the distributions of local charges and 
magnetic moments arising from inhomogeneity in the site environments. 

Even at the mean field level of UHF, the combined effects of disorder and electron 
interactions in leading to strongly inhomogeneous ground states yield rich behaviour, as 
shown by a recent study of the zero-temperature UHF phase diagram of a site-disordered 
Anderson-Hubbard model0 (AHM). Further, at RPA level, these ground states in turn give 
rise to strongly inhomogeneous magnetic response properties!^!. In particular, disorder is 
found to lead to a significant enhancement in local static susceptibilities over those in the 
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pure system, this enhancement being strongly site-differential. 

The behaviour of static and dynamic susceptibilities in both the spin and charge channels 
ultimately reflects the character and distribution of collective excitations about the ground 
state. The observed site-differential enhancements mentioned above naturally raise questions 
both of how disorder affects the excitation spectra of the model and of the spatial distribution 
and localization characteristics of the excitations. In addition, the stability of the long- 
ranged order present in the mean-field ground state towards zero point spin fluctuations, 
and the low-temperature properties of the model, are also governed by the RPA excitations. 

We here investigate of the RPA spectrum of a d — 3 Anderson-Hubbard model with 
Gaussian site-disorder, throughout the disorder-interaction phase plane. In Section II, we 
discuss generic properties of the model and describe the procedure for calculating and char- 
acterizing the collective excitation spectrum at RPA level. After a brief discussion of the 
pure Hubbard model in Section III, we focus in Section IV on the AHM, concentrating on 
the distribution, character and localization properties of the low-frequency excitations, with 
particular emphasis on the existence or otherwise of a spin wave-like band in the transverse 
excitation spectrum at finite disorder. We conclude in Section V with a discussion of the 
extent to which a recently developed mapping of the excitations of the Hubbard mo 
onto those of an effective spin model, can be applied to the disordered system. 



II. COLLECTIVE EXCITATION SPECTRUM AT RPA LEVEL 

The Hamiltonian we consider is given by 

H = ein ia - t c \<r c i<r + W n iv n i-v (2-1) 

where t is the hopping matrix element, U is the (repulsive) on-site Coulomb interaction and 
the (ij) sum is over nearest neighbour sites on a d = 3 simple cubic lattice. The site energies 
{ej} are drawn randomly from a common Gaussian distribution g(e) of variance A 2 which, 
together with Eq. Q2.1|) , specifies the Anderson- Hubbard model considered here; and we 
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focus exclusively on half-filling. Making a spin-rotationally invariant UHF approximation to 
the interaction term in Eq. fl2.1|) allows the Hamiltonian to be expressed as H = H + Hi 
with 

i,a (ij),a 

+Uj2{hn i n i -2S i -S i } , (2.2) 

i 

where Si is a spin-| operator and the overbar denotes an expectation over the UHF ground 
state. The zero-temperature phase diagram of H in the (A/t,U/t) plane has been discussed 
extensively in Ref. [11]; aspects of this study necessary for the present work will be summa- 
rized briefly in §IVA. Here we note that, while a rich variety of phases arises, all mean-field 
magnetic ground states are found to be Ising-like (i.e. fully collinear local moments, with 
Si X — — S iy for all sites i), with 5'* ot = Y^i Siz — 0. Collective excitations about these 
ground states, obtained via the RPA, then decouple into two setsB for frequencies cu > 
there exist N 2 /2 spin excitations transverse to the local moment z-axis, and N 2 /2 longitu- 
dinal spin and charge excitations. In the non-disordered (A/t = 0) Hubbard model there 
is a gap to longitudinal spin excitations for all U/^, reflecting the presence of a gap in 
the single-particle density of states (DoS). By contrast the transverse excitation spectrum 
is gap less (as reflected by the presence of Goldstone modes). With disorder present, the 
dominant low-energy collective excitations are again found in the transverse spin channel. 
We thus focus here on the transverse spin spectrum, since it is these excitations which 
govern measured susceptibilities (the existence of a spin-flop transition implying that only 
transverse excitations are in practice probedlil), and which are dominant in determining 
low-temperature thermodynamic properties. 

A. Determination of the transverse excitation spectrum 



The equations of motion for the RPA particle-hole excitations are straightforwardly de- 
rived (see e.g. Ref. [L6|). However, solution of these equations requires diagonalization of a 
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2N 2 x 2N 2 non-Hermitian matrix, a procedure limited to impractically small system sizes. 
An alternative — and viable — procedure is to exploit the fact that, because the interaction 
term H\ contains only on-site terms, expressions for the RPA dynamic susceptibilities (de- 
rived either diagrammatically or via linear response theory^), involve only matrices of order 
N. Specifically, the transverse spin susceptibility 



% J dt e^(0\T{Sr (t)S f}\0) RPA , (2.3) 



is given by 



X - »= V» l-t/ u %- + H ■ , (2.4) 



-i 



where the N x N matrix °x + ( UJ ) is the corresponding UHF transverse susceptibility ob- 
tained by replacing |0)rpa by |0)uhf m Eq. ( |2.3| ). This is in turn given explicitly by 



j dia1 a ja] a i0l a j/3l 



a>F>p I E ^ -E^-u-vq 



E al - + u - if] J ' 

Here F is the Fermi level, 77 = 0+ and the {E aa } are the UHF single-particle energies, 
with eigenvectors \^ aa ) = J^i^aal^ia) expanded in a site basis (and which are pure spin- 
orbitals for Ising-like UHF solutions^). We add in passing that the Fermi energy is pinned 
at Ep = ~U for all disorder strengths A > 0, since the familiar particle- hole symmetry 
characteristic of the non-disordered limit is preserved for all A > with the symmetric 
Gaussian g(e )§. 

Eq. ( |2.5| ) is just the Lehmann representation for the retarded polarization propagator 
appropriate to the UHF ground state, with poles at the Stoner excitation energies. When uj 
does not coincide with a pole, this real symmetric matrix is diagonalized by an orthogonal 
matrix V(u;) with eigenvalues {A 7 (u;)}. Most importantly, as is evident from Eq. fl2.4| ), this 
transformation also diagonalizes X~ + ( UJ ) with corresponding eigenvalues A 7 (u;)/[1 — UXy(u))], 
viz: 
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RPA transverse collective excitations correspond to the cu-poles of x h ( Co ')) which thus occur 
at 1 — UXyiu) = 0, i.e. whenever an eigenvalue of °x l ~(^) crosses 1/U. Further, since 
xffiyj) = xJi + (~ u>), a knowledge of the (positive and negative) poles of x~ + ( UJ ) alone is 
sufficient to determine the full transverse spin excitation spectrum. 

Two final points concerning the poles X~ + i UJ ) should be made before proceeding. First, 
the broken symmetry of any magnetic UHF state implies the presence of a zero-frequency 
Goldstone mode in both x~ + an d X + ~-> ensuring that the transverse excitation spectrum 
is gapless. Second, the tu-poles occur correctly on the real axis only if the UHF ground 
state is fully stable against particle- hole excitations (see e.g. Ref. [17|). An eigenstate of 
Hq which is not a true minimum on the UHF total energy surface leads to imaginary-^ 
excitations which carry the system away from the saddle point, rendering it unstable. In 
this regard, and even when the mean-field solutions are Ising-like, use of a fully unrestricted 
and spin-rotationally invariant Hartree-Fock approximation is an essential prerequisite, since 
it ensures the mean-field states are minima. 

In studying collective excitations at RPA level, a pole search on x~ + ( UJ ) ( Eqs (2.4,5)) 
is powerful, since it allows study of much larger system sizes than direct solution of the 
RPA equations. We now outline this procedure. From Eq. (|2.5|) , °x~ + ( UJ ) itself has poles 
at the UHF excitation energies but, in the intervals between the poles, the eigenvalues 
{A 7 (u;)} of °X~ + ( UJ ) are continuous functions of uj. For a given disorder realization, we thus 
determine the excitation spectrum as follows, (i) Solution of the UHF problem yields the 
UHF excitation energies and hence the intervals over which the eigenvalues A 7 (u;) are finite, 
(ii) Diagonalization of °X~ + { UJ ) f° r u a ^ t ne beginning and the end of each interval yields 
the number of eigenvalues which have crossed 1/U — and hence the number of RPA poles — 
in each interval, (iii) Finer resolution may then be obtained using, for example, a bisection 
technique. 
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B. Character of the transverse excitations 



Important information concerning the character of a given excitation of frequency uo p can 
be obtained from the eigenvector coefficients {V^iup)}. 

In the strong coupling limit U/t oo, the resultant AF Heisenberg model contains 
N transverse particle-hole excitations, corresponding to linear combinations of solely on- 
site spin flips. These are just the familiar pure spin waves. At finite U/t, off-site spin 
flip excitations also contribute to the transverse spectrum, such that while there is spin- 
charge separation (in the sense that the transverse spin and longitudinal /charge channels 
decouple^) , the transverse excitations themselves possess some charge-transfer character. 
It is particularly instructive to evaluate the on-site spin flip (or spin wave) contribution to 
a given excitation, as now considered. 

The Lehmann representation for the diagonal elements of the RPA susceptibility is given 
generally by 

n£nO-W-iri £ n0 + LJ- IT] 

where |0) and \n) are RPA ground and excited states and £ n o are the RPA excitation 
energies. From Eq. ( |2.6| ), Xu + ( UJ ) is expressible in terms of the eigenvalues and eigenvectors 



of° X -+(u<): 



*;>) = XXMr^) (2.8) 



To compare Eqs Q2.7|) and (|2.8| ), we analytically continue Eq. ( p.8|) by the replacement 
A a (a>) — > \ a (u> + is), where s is an infinitesimal. Consider any particular pole, u p , of Eq. 
(p.8|) , such that \{oj p ) — 1/U. For frequencies u ~ u p , is dominated by the a = 7 

term and is given to leading order by 

Xi +{u~u p )= "> . (2.9) 

Comparison of Eqs Q2.9| ) and (|2.7|) implies 



sgn (^) = sgn(cj p ) , (2.10a) 

\ / Up 

and 

sgn(s) = sgn(^p) (2.10b) 

(in order that the poles of x^ + ( UJ ) occur in the correct half-plane); together with the central 
result that, for the given pole, 
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K0|ff|n)l'= V&u,) (2.11) 



u 2 



where p = — sgn(a; p ). Further, since the eigenvector coefficients {Vi^(u))} are by construct 
normalized, the sum rule 

E\m\n)\ 2 = — 4— (2-12) 



U 2 



du> 

follows. That this sum is not unity simply reflects the fact that, at finite U/t, RPA transverse 
spin excitations have finite weight outside the on-site spin flip subspace, i.e. matrix elements 
of the form (0|cjiCjf|n) with i ^ j are non-zero; for example it may be shown that the general 
RPA normalization expressed in a site basis implies that J2ij I (0| c !| c jTl n ) | 2 = 1 f° r RPA 
excitations corresponding to positive-u; poles. Only in the strong coupling limit U/t — > oo 
are these off-site excitations effectively projected out of the transverse spin spectrum. 

The interpretation of the eigenvectors {V 7 (u;)} is therefore clear: from Eqs ( p.llj ) and 
Q2.12j ) the set of coefficients {V^iup)} corresponding to a particular A 7 (co> p ) = jj describe the 
spatial distribution of the on-site spin flip component of the collective excitation, while the 
prefactor VL{oj p ) = l/(U 2 \dXy/du\ u ) gives the weight of the RPA excitation in the subspace 
of on-site spin flips. Excitations with Q(u) p ) ~ 1 are predominantly spin waves, while those 
with Q(uj p ) 1 have significant off-site 'charge transfer' character. (Note that, to generalize 
to a disordered system, we use the term 'spin wave' to denote any excitation with Q ~ 1, 
regardless of its ^-resolution.) 
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The above interpretation may also be exploited to define an inverse participation ratio 
(IPR) for the collective excitations, in analogy to the familiar IPR employed in the study of 
localization of single-particle states^, via 

L(u p ) = J2V*(u p ) . (2.13) 

i 

Using Eq. (2.12), this may be expressed as 

L M = 7^ | /n , QP | \, 2 x2 :p = -sgn(u; p ) . (2.14) 

L(u) p ) is on the order of the inverse of the number of sites overlapped by the excitation 
at uj p : for an excitation receiving uniform contributions from m sites, L(u) ~ 1/m. In 
the thermodynamic limit, L(u) is thus zero for a delocalized excitation and non-zero for a 
localized excitation. For a finite-sized system L > necessarily, requiring the estimation of 
a threshold IPR, appropriate to a given system size. While for the single-particle excitations 
of the (non-interacting) Anderson model, such a threshold IPR can be addressed via finite- 
size scaling^, we have no information on the A-scaling of L(u) in the present system. In 
§IV, however, we argue that a qualitative understanding of the localization characteristics 
of transverse spin excitations may nonetheless be obtained. 



III. EXCITATION SPECTRUM OF THE PURE HUBBARD MODEL 

Before discussing the effects of disorder on the transverse excitation spectrum of the 
Anderson-Hubbard model, we first consider briefly the important non-disordered limit. In 
this case, for all U/t > 0, the UHF single-particle spectrum D(E) = \ J2aa ~ E aa ) has a 
band gap of magnitude U\fi\, with \fi\ the self-consistent UHF local moment magnitude; and 
in which lies the Fermi level Ep = Consequently, the pure UHF particle- hole excitation 

spectrum consists of a single Stoner band beginning at u = U\fi\, with a maximum at u ~ U. 
The full RPA spectrum of collective transverse spin excitations consists, by contrast, of two 
distinct components (we consider u > throughout, since the spectrum is symmetric in uj). 
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First, a low frequency band containing precisely N excitations of predominantly spin wave 
character (discussed below), and extending down to uj = as implied by the presence of the 
Goldstone modes for all U/t > 0. Second, a high frequency band containing (N/2)(N — 2) 
Stoner-like excitations, which closely resembles the corresponding pure UHF transverse spin 
spectrum: it begins at uj = U\[/,\ — the gap in the UHF single-particle spectrum — and has 
maximum spectral density for u ~ U. This Stoner-like band is effectively eliminated as 
U — > oo. For all U/t > 0, however, it may be shown to be separated from the top of the low 
frequency band by a spectral gap, i.e. there is a persistent separation of scales between the 
low energy spin wave-like band and higher energy Stoner-like excitations. For this reason 
we focus exclusively on the low-o; component of the spectrum. 

The prescription of §11, reformulated in a two-sublattice basisB appropriate to the half- 
filled, non-disordered model, allows the RPA spectrum to be determined; in practice we 
find that a lattice of 32 3 q-points (corresponding to N = 32768) is required to give a 
sufficiently dense spectrum. Fig. [I] shows the resultant density of the lowest iV transverse 
spin excitations, given by 

N t (uj) = N- 1 5 ( u ~ °°p) 

a.x 

where co> max is the frequency of the N th excitation, for large coupling U/t = 20. Also shown 
is the corresponding spectrum for the AF Heisenberg model with solely nearest-neighbour 
(NN) J = 4t 2 /U, given by N^ sw (uj) = N' 1 £ q c^ - c^ q ) with! 

^ = ±37^1-72 (3.2) 

where 7 q = | ELi cos(q a a) and a is the lattice constant. The spectrum of low- frequency 
transverse spin excitations of the Hubbard model thus closely resembles that of the NN 
Heisenberg model at large U / 1 — as expected, since as U/t —>■ 00 the former maps exactly 
onto the latter — although as seen there is a not inappreciable quantitative difference even 
for U/t = 20. 

Less predictably, this qualitative similarity persists down to weak coupling. However, as 
U/t is decreased, the energy scale for low-frequency RPA excitations rapidly deviates from 
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that of the NN Heisenberg model with J = 4t 2 /U. Fig. |2] shows the effects of gradually 
decreasing U/t down to U/t = 2. Initially, the spin- wave band is pushed to higher energy, 
consistent with the strengthening of NN exchange couplings and reflecting also the increasing 
importance of beyond-nearest neighbour couplings as shown in a recent paper0. Below 
U/t ~ 5, however, this band starts to migrate to lower energies, reflecting an increasing 
repulsion with the approaching Stoner-like band. 

To investigate quantitatively the spin-wave character of the excitations, we calculate 
fl(uj p ) = 1/ (U 2 \d\/ du\ Wp ) for uj p < u; max . The lowest excitation is of course the Goldstone 
mode corresponding to a global spin rotation, and is a pure spin wave for all U/t > 0: Q(uj p = 
0) = 1. Since excitations at the upper edge of the spin wave-like band are naturally most 
strongly influenced by the encroaching Stoner-like band, it is found that D,(u p ) > Q(uj mauX ) 
for all Up < u max . The first column of Table I thus shows the [7/t-dependence of Q(u; max ). 
From this it is seen that for U/t > 3 the excitations are predominantly spin waves - in fact 
for U/t > 5 they are in essence pure spin waves, Q = 1 — commensurate with a significant 
gap between the spin wave and Stoner-like bands. For U/t < 2 — 3, this gap diminishes 
rapidly and there is a substantial loss of spin wave character for excitations with uo ~ uj max . 
Nonetheless, as is evident from the above, the low energy transverse spin excitations of the 
Hubbard model remain predominantly spin waves for all but the lowest U/t. 



IV. ANDERSON-HUBBARD MODEL 

Before considering the effects of disorder on the RPA spectrum, we first briefly outline 
relevant properties of the inhomogeneous UHF ground state, since it is particle-hole excita- 
tions out of this state which form the basis for the RPA collective modes; a full account of 



the mean- field phase diagram is given in Ref. 11 
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A. Mean-field ground state 



On the introduction of disorder, UHF solutions remain Ising-like, with Sl ot = 0. Three 
distinct magnetic phases are found: disordered paramagnetic (P), antiferromagnetic (AF) 
and spin glass-like (SG) phases. Both metallic (M) and insulating (I) phases also occur, with 
the dominant metal-insulator transition (MIT) M— >gapless I, driven by Anderson localiza- 
tion of single-particle states at the Fermi energy Ep — |C/. We now consider two specific 
paths across the disorder-interaction phase plane, which will be employed in the next section 
to illustrate the combined effects of A and U on the RPA transverse spin spectrum. 

First, consider fixing U/t = 12 and varying disorder, A/t. As A/t is increased from 
zero, disorder leads to broadening and eventual overlap of the Hubbard bands in the UHF 
single-particle spectrum, such that Ep lies within a pseudogap in the total DoS D(E). 
Studies of the IPR L'(E) = J2i \ a iaa\ 4 for single-particle states of energy E aa = E, with 
E ~ -Ef^ imply that the system passes directly from a gapped to gapless insulator in 
this region of the (A/t, U/t) phase plane (at A/t ~ 2.5). Fig 0(A) shows D(E) and the 
IPR profile L'(E) averaged over N = 512 site disorder realizations for [A/t, U/t) = (3, 12). 
A pronounced pseudogap is evident in D(E), in which Ep lies. Correspondingly, L'(E) 
shows a maximum at E — Ep and minima for states in the centre of the Hubbard bands 
where D(E) is maximal. Using a threshold IPR appropriate to the chosen system size0 
shows that Fermi level states are localized: the system is a gapless Mott- Anderson insulator 
(and remains thus with increasing A/t at the chosen U/t = 12). Analysis of the Fourier 
transform of the local moment S z (k) = N" 1 ^^ shows however that the system 
remains antiferromagnetically ordered: although increasingly 'dirtied' by disorder, |5 z (7r)| ~ 
\fj,\ (where the mean moment per site — N~ l Y^i l/^l); with little weight at other k- vectors. 
Note however that while the phases of the moments remain locked in AF alignment, the 
distribution of local moment magnitudes over the sites is highly disordered. This is evident 
from the site-resolved distribution of local moment magnitudes, |/x(e)| = iV" 1 Y J i-.e i =e | A*» I j 
with N e the number of sites with e« = e; as shown in Fig 0(B), again for (A/t, U/t) = (3, 12). 
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Sites with energies |e$| < \U, disposed randomly on the lattice, carry strong atomic-like 
moments, while sites with lower or higher site energies have \m\ ~ 0, being respectively 
largely doubly occupied by electrons or empty. This profile is found to vary little with 
disorder for fixed U/t = 12, aside from minor erosion of local moment magnitudes: the 
primary effect of A/t is to determine, via the site-energy distribution g(e), the fraction of 
sites (around e = 0) which carry strong local moments. 

Also shown in Fig. 0(B) is a measure of the distribution of quasiparticle states of given 
energy E over the sites: H(e; E), given by0 

H ^ E >- D(E) ' (41a) 

where 

D(e; E) = £ £ \a ia(7 \ 2 5(E - E aa ) (4.1b) 

i:ei=e out 

is the partial density of states for sites of bare site energy e. H(e; E F ) is thus the quantum 
probability density that electrons in Fermi-level quasiparticle states will be found on sites 
of bare site energy e. As discussed in Ref. [II], and exemplified in Fig. 0(B), H(e;E F ) is 
peaked strongly at the boundaries of the local moment distribution |e| ~ \U] i.e. sites with 
strong local moments participate but weakly in Fermi-level quasiparticle states, which are 
in contrast dominated by sites with bare site energies close to the local moment boundaries. 
Sites with strong local moments, around e = 0, are in turn found to contribute mainly to 
quasiparticle states with E ~ and E ~ U, i.e. around the maxima in the single-particle 
spectrum. 

The second path considered is fixing the disorder at A/t = 3, and decreasing the scaled 
interaction strength U/t from U/t = 12. Note that, since the U/t — > oo limit of the disor- 
dered Hubbard model for any A/t is the non- disordered AF Heisenberg model, decreasing 
U/t may be thought of as effectively increasing (by 'switching on') the role of disorder. As 
U/t is reduced, the fraction of sites at the local moment boundaries ~ g(±^U) — which 
give the dominant contribution to Fermi-level states (Fig. 0(B)) — increases. The system 
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remains 



is ultimately driven metallic at a critical U/t ~ 7. 20. However, magnetic ordering; 
characterized by a strong peak in S z (k) at k = n until U/t ~ 2.5, whereupon a sharpll! 
transition occurs to a spin glass-like state, with S z (k) receiving essentially uniform 0(1/N) 
contributions from many k vectors; see Ref. [12]. In this phase, although the mean local mo- 
ment per site is typically small (~ 0.05), relatively strong local moments are in fact found 
on a small fraction of sites with |e| < ~U for a given disorder realization. As discussed in 



Ref. [11], this reflects the disorder-induced production of statistically rare local environments 
in which atomic-like moments are not significantly eroded by electron hopping processes. 



B. Collective excitation spectrum of the Anderson-Hubbard model 

We now turn to the effects of disorder on the transverse RPA excitation spectrum. The 
principal questions we seek to address are: 

• What is the effect of disorder on the density of transverse spin excitations? Does a 
low-energy spin wave-like band exist for finite disorder, and what is the character of 
the associated excitations? 

• What is the effect of disorder on the localization characteristics of low-lying collective 
excitations? 

• What is the distribution of transverse spin excitations of given energy over the sites, 
and how does this relate to the inhomogeneous distribution of local moments in the 
mean-field ground state? 

With disorder present, there are no simplifying symmetries which allow very large system 
sizes to be studied, as in the A = case. However, this lack of symmetry — specifically, the 
absence of degeneracies in the UHF spectrum — enables reasonable statistics to be gained by 
sampling many disorder realizations for system sizes much smaller than those required for 
A = 0. Numerical evaluation of the entire transverse excitation spectrum — including the 
Stoner-like excitations — remains nonetheless a large problem: a full pole search for a system 
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of N sites requires diagonalization of a. N x N matrix in the vicinity of |iV 2 poles, and is 
therefore 0(N 4 ) for each disorder realization; an ensemble of which needs subsequently to be 
sampled. For calculations of the full transverse spin spectrum we have thus chosen N = 64 
site systems, which are found to be sufficient for reproduction of its general features; while 
for more detailed study of the low energy portion N = 216 is employed. 

For fixed disorder A/t = 3 and various values of the interaction strength U/t < 12, Fig. 
§ shows the full RPA density of transverse spin excitations, together with the corresponding 
UHF particle-hole excitation spectrum (the 'pure Stoner' spectrum). At U/t = 12, where 
the mean-field ground state is a gapless AF insulator, a pronounced spin wave-like band is 
again evident in the RPA spectrum at low lu, as for the non-disordered limit (although the 
pure Stoner spectrum now extends down to uo = 0); the remainder of the RPA spectrum 
closely resembles the Stoner spectrum. On decreasing U/t, the pure Stoner spectrum shifts 
to lower energy, since the separation between the maxima in the corresponding single-particle 
spectrum D(E) decreases and the pseudogap at E F is gradually eliminated. As a result, 
the spin wave-like band is gradually shifted to lower energy ('softened') by, and ultimately 
absorbed into, the Stoner-like band. By U/t = 3, close to the AF-SG border, the pure 
Stoner and RPA spectra appear virtually indistinguishable. 

On fixing U/t and increasing disorder, a somewhat different picture emerges. Fig. || 
shows the RPA and pure Stoner spectra for U/t = 12 and various values of A/t. At 
low disorder A/t = 1, for which the UHF ground state is a AF Mott-Hubbard (gapped) 
insulator, a spin wave-like band, well separated from the remainder of the excitations, is 
clearly distinguishable. Closer resolution shows that this band is in fact shifted to slightly 
higher frequency compared to the non-disordered case, despite the intuitive expectation that 
disorder should lead to a 'softening' of the spin waves. This upward shift is correspondingly 
manifest in an initial decrease in the uniform static susceptibility \u = J2ij Xij + (0) upon 
the introduction of disorder: as follows from the Lehmann representation of the susceptibility 
Eq. (|2.71), a decrease in density of low- frequency excitations with significant weight in the 



on-site spin flip subspace leads to a decrease in Xu- On increasing disorder further, the UHF 
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single-particle gap closes and the pure Stoner spectrum acquires weight down to 00 = 0, Fig. 
[5]. The maximum in the Stoner-like band in the RPA spectrum — occurring at uj ~ U— 
changes little with disorder, but its disorder-induced broadening increasingly repels the spin 
wave-like band which is thus by contrast shifted to lower frequency (see Fig. |^). This is 
the 'softening' of spin waves with disorder, which in turn is reflected in an increase in the 
uniform static susceptibility Xu- Fig. || shows the disorder- averaged Xu vs A/£ for U/t = 12 
(and we add that in the non-disordered limit, the resultant Xu for the N = 64 site system 
reproduces accurately the N oc value, which can be obtained analytically). The above 
behaviour is clearly demonstrated: after an initial slight decrease, Xu progressively increases 
with increasing disorder as the spin wave-like band moves to lower frequencies. Note further 
that increasing disorder for fixed U/t has a considerably less pronounced effect on the spin 
wave-like band than decreasing U / 1 for fixed disorder: decreasing interaction strength shifts 
the maximum in the Stoner-like band to lower energy (Fig. while increasing A/t by 
contrast merely broadens the band, as shown in Fig. |5]. 

The effect of disorder on the position of the spin wave-like band also provides a natural 
explanation for the strongly site-differential enhancements of local susceptibilities with dis- 
order, as found in a recent RPA study of the present mo deli. Specifically, the total unform 
susceptibility is deconvoluted via 

Xu = J deg(e) X (e) (4.2) 

where x( e ) — ^7 1 Y.i:e i =eXi with Xi — XwXy(O). x( e ) is thus the mean local susceptibility 
for sites of given site energy e. In ref. |13] (Fig. 2) it was found that an increase in disorder 
leads to a strong enhancement of x( e ) for sites in the local moment range |e| < ^U, while 
sites outside this range, whose susceptibilities have typical Pauli-like values, are by contrast 
only weakly affected by disorder. A site-differential analysis of the RPA excitations forming 
the spin wave-like band, in analogy to H(e;E) for single-particle states (Eq. ( |4.1| )), shows 
that these excitations have weight almost exclusively on the strong local moment-carrying 
sites, and receive little contribution from sites outside the local moment range (consistent 
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with their interpretation as 'spin wave'-like). As discussed above, the principal effect of 
increasing disorder is ultimately to shift the spin wave-like band to low frequency, while the 
Stoner-like band is merely broadened. The disorder- induced 'softening' of the spin wave-like 
band therefore manifests itself almost exclusively in x(e) for strong moment-carrying sites, 
as found in Ref. O; and has little effect on susceptibilities outside the local moment range, 
which are instead dominated by Stoner-like excitations. 



C. Character of the transverse excitations 

From the preceding section, it is clear (Figs |],[5]) that a low-energy feature in the density of 
RPA transverse spin excitations indeed exists in the AF region of the (A/t, U/t) phase plane. 
However, the preceding analysis by itself gives no information on the spin wave-like character 
of the excitations which form this band, or with which sites they are associated. To probe the 
nature of these low-energy excitations, we have performed pole searches using N = 216 sites 
for many disorder realizations at a selection of points in the (A/t,U/t) plane, focusing (as 
in the non-disordered limit) on the lowest N excitations, with fl(a> p ) = l/(^ 2 |<9A/<9u;| w=Wp ) 
evaluated at each RPA pole. For illustration, the second and third columns respectively of 
Table I show Q(u) averaged over the lowest N RPA excitations for (a) fixed A/t = 3 and 
various interaction strengths U/t; and (b) for fixed U/t = 12 as a function of disorder. 

For fixed U/t = 12, it is seen that while increasing disorder leads to a diminution of spin 
wave character, the effect is relatively weak and the low energy transverse spin excitations 
retain strong spin wave character. This is consistent with the fact that the system remains an 
AF insulator, and with the persistence (Fig. 5) of the pronounced low energy spin wave-like 
band in the transverse spin spectrum. By contrast, reducing U/t at fixed disorder leads to a 
more pronounced reduction in spin wave character, consistent with the progressive erosion 
of the spin wave-like component in the RPA transverse spin spectrum (Fig. 4) discussed 
above. Again, this is physically natural: with decreasing interaction strength the system 
undergoes a (gapless) insulator— >metal transition at U/t ~ 7.2 and an increase in off-site 
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character in the low-energy spin excitations is thus to be expected. Note however that, even 
at low U/t, the RPA excitations have substantially more spin wave character than pure 
Stoner excitations in the same frequency range, as may be verified by considering the weight 
of the latter in the on-site spin flip subspace, given by 

OhpW) = E E *L,*%-.&H - [E„ - E^]) (4.3) 

icr a>f3 

This is found to be an order of magnitude smaller than the corresponding for low- 

frequency RPA excitations. 



D. Localization characteristics of 1ow-cj excitations 

Loss of translational symmetry on the introduction of disorder is expected to lead to 
localization of some or all collective excitations, in analogy to Anderson localization of 
single-particle states. Consider the IPR L(u p ) of a collective transverse spin excitation 



of energy u p , as defined in Eq. ( 2.13Q . To distinguish between localized and delocalized 
excitations a threshold IPR appropriate to a given system size should be established, such 
that excitations below (above) this value are deemed extended (localized). In principle, this 
could be obtained by finite-size scaling, but in practice we are limited to system sizes on the 
order of N = 216. For single particle states of the non-interacting Anderson model, finite- 
size scaling yields a threshold IPR L' c ~ 0.1 for N = 216 sites. We correspondingly adopt 
a threshold L c ~ 0.1 for collective excitations of the Anderson-Hubbard model. This value 
itself should not of course be taken very seriously; it simply serves as a useful qualitative 
guide in the following discussion. 

Fig. [7| shows the IPR L(u) of the lowest energy excitations for U/t = 12 and various 
values of the disorder. For weak A/t = 1.5 (Fig. [7] A), the majority of excitations are 
delocalized, as expected, with a localized 'tail' at high u. Increasing disorder to A/t = 
3 and 4.5 leads to the occurrence of further sparse, localized excitations at very low uj, 
superimposed over a background of delocalized excitations (Fig. [7] B,C); and the density of 
localized excitations rises on further increasing disorder (Fig. [7| D). 
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The above behaviour is clarified by considering the corresponding IPR probability dis- 
tributions, P(L) (Fig. |8p. This shows a peak at low L, corresponding to the 'delocalized' 
backbone, with a long tail to high L, reflecting the presence of localized excitations. For 
low A/t = 1.5, the peak is sharp and centred well within our approximate criterion for de- 
localized excitations. As A/t is increased, the distribution broadens and the most probable 
value of L shifts to higher values. 

To understand the above behaviour, we site-resolve excitations of given IPR L, in analogy 
to the site distribution of quasiparticle states, if (e; E) (Eq. (|4.1|) ), via the quantity 

%AM) -nur (4 - 4) 

(such that / deH RP j±(e; L) = 1). Here, N RPA (L) = iV -1 £(£ — L(u p )) is the density of 
excitations of fixed IPR L, and iV R p A (e; L) is the partial density of such excitations on sites 
of given site energy e, given by 

iV RPA (e;L) = iV ( r 1 £ £ \V^(u p )\ 2 6(L - L{u p )). (4.5) 

Hnp\(e; L) then gives the overlap of excitations of given IPR L on sites of different site 
energies. For (A/t, U/t) = (3, 12), Fig. |9] shows the e-dependence of H RPA _ averaged over all 
L < 0.1 and L > 0.3, to illustrate the site-differential character of extended and localized 
excitations respectively. Localized excitations are seen to be dominated by sites close to 
the local moment boundaries |e| ~ U/2, while sites in the range |e| < U/2, which carry 
large moments, make the dominant contribution to delocalized excitations. As expected, 
sites outside the local moment regime, which possess very small moments, do not contribute 
significantly to any of the low energy excitations; rather, they participate in Stoner-like 
excitations of significantly higher energy. 

The origin of the localized and delocalized low energy transverse spin excitations is then 
readily rationalized. Sites with site energies in the range |e| < U/2, which possess large local 
moments, contribute most significantly to single particle states within the Hubbard sub- 
bands. Since these states are the most delocalized, collective excitations involving them will 
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likewise be delocalized. At UHF ('pure Stoner') level, transverse spin excitations between 
single-particle states within the Hubbard bands incur an energy cost of order U ; but in the 
RPA spectrum, via inclusion of particle-hole correlations, this energy is renormalized down 
to a spin wave-like scale on the order of an exchange coupling. This is the origin of the 
delocalized low-energy transverse spin excitations. By contrast, sites with |e| ~ U/2 close 
to the local moment boundaries participate dominantly in single-particle states close to the 
Fermi level, or an energy ~ U above or below it. For U/t = 12 such states are strongly 
localized, hence so too are collective excitations involving them. Again, such excitations are 
renormalized to a low-energy scale via the particle-hole interactions inherent in the RPA. 
Both localized and delocalized transverse spin excitations thus arise at low energies, reflecting 
the underlying localization characteristics of single particle states over a very wide energy 
range in the associated single-particle spectrum. Note further that these arguments do not 
imply coexistence of localized and extended transverse spin excitations at the same energy 
in the spin wave-like band; and the apparent occurrence of such in eg Fig. 7D is largely 
a consequence of sampling different disorder realizations. While the present, qualitative 
considerations do not of course permit identification of mobility edges in the spin wave-like 
band, it is at least plausible from Fig. 7 that such exist. 

From the above, the system can thus be regarded as composed of two essential com- 
ponents: an AF large-moment component which gives rise to the backbone of delocalized 
excitations, and a set of sites carrying smaller moments which give rise to localized excita- 
tions. As disorder is increased, the fraction of sites in the former category decreases, leading 
to a shift in the peak in P(L) (Fig. |8]) to higher L and a concomitant broadening of the 
distribution. The probability of sites with |e| ~ U/2 correspondingly increases, and hence 
localized excitations ultimately become the dominant theme. 

The above discussion is confined to a region of the (A/t,U/t) phase diagram in which 
the mean field ground state is an AF insulator. What of the metallic AF and SG phases? 
Fixing A/t = 3 and decreasing U/t progressively from U/t = 12 down to U/t = 6 leads 
(at U/t ~ 7.2) to a transition from an AF insulator to an AF metal. However, the IPR 
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profile of the low-energy transverse spin excitations shows essentially no qualitative change 
from that for U/t = 12: a backbone of delocalized excitations remains present, with local- 
ized excitations superimposed at low and higher u/t. In the metallic region, however, the 
localized excitations now arise principally from transfer of a spin between a delocalized UHF 
state near the Fermi level to a localized UHF state ~ U above or below Ep. The localized 
state overlaps only a few sites, while the delocalized state has weight on a finite fraction 
of sites. The excitation thus appears localized in the IPR profile since the on-site spin flip 
contribution (which is of course the component of the excitation probed, see Eq. (2.14)) can 
occur only from sites overlapped by the localized state. 

On decreasing U/t further to U/t ~ 2.5, the mean field ground state becomes a SG 
metal. Here, the RPA and pure Stoner excitation spectra appear virtually identical (see e.g. 
Fig. |j). It is indeed found that the majority of low energy RPA excitations are little more 
than weakly renormalized Stoner excitations, as expected from the low value of Q in Table 
I. However, at very low energies (cu < 0.3t), this is not the case. As in the AF metal regime, 
low energy localized modes appear in the RPA spectrum, while the pure Stoner spectrum at 
such low energies is composed of entirely delocalized modes (since the system is metallic). 
These localized excitations are found to have dominant weight on those rare sites in the spin 
glass that possess substantial local moments, flips of which would incur an energy cost of 
~ U within UHF, but occur at significantly lower energies in the RPA spectrum. 



V. DISCUSSION 

In recent papersB'El, we have described an approximate mapping of the low-frequency 
transverse spin excitations of the non-disordered Hubbard model at finite U/t, on to those 
of an effective underlying Heisenberg model. The [7/t-dependent effective exchange cou- 
plings are not restricted to purely nearest neighbour interactions; they follow solely from a 
knowledge of the RPA excitations about the broken-symmetry UHF state at zero temper- 
ature. The mapping becomes exact as U/t — > oo and, conjoined with an Onsager reaction 
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field approach, enables extraction of thermodynamic properties in the thermal paramagnetic 
phase. The resultant Neel temperatures, spin correlation functions and magnetic suscep- 
tibilities were0E3 found to be in good agreement with Quantum Monte Carlo calculations 
over a wide U / t-range from strong to weak coupling. We now assess the extent to which the 
present findings support application of such a mapping to the Anderson-Hubbard model. 



As discussed in Refs |1J, [15], the mapping is accurate provided there is a persistent separa- 
tion between low energy spin wave-like excitations and higher energy Stoner-like excitations; 
i.e. provided in practice a discernible spin wave-like band, comprising excitations with an 
appreciable degree of on-site spin flip character, is present in the full RPA transverse spin 
spectrum. As shown in §111, this indeed holds for the pure Hubbard model for all U/t > 2 — 3. 
For the Anderson-Hubbard model, Fig. |5| and Table I show that this also holds for all dis- 
order strengths studied, at relatively large U/t = 12. Thus, the mapping should certainly 
be accurate well within the disordered AF phase of the model. On fixing disorder and re- 
ducing interaction strength towards the AF-SG boundary, the spin wave character of the 
low-frequency band decreases significantly and the band itself is ultimately absorbed by the 
Stoner-like portion, such that for A/t = 3 the mapping would not be quantitatively accurate 
for U/t ~ 4 and below. Not unexpectedly, therefore, the mapping would not be applicable 
to the low-U/t SG phase; but we do anticipate it to be accurate throughout the major part 
of the disordered AF phase in the (A/t, U/t) plane, wherein the effect of disorder on, for 
example, the Neel temperature can be investigated. This is the subject of a forthcoming 
papeiS. 
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FIGURES 

FIG. 1. Density of low- frequency RPA excitations Nt(w) for the pure Hubbard model with 
U/t = 20 (solid line), compared to the linear spin wave spectrum of the pure AF Heisenberg model 
with nearest-neighbour J = 4t 2 /U (dashed line). 

FIG. 2. Density of the N lowest-frequency RPA excitations Nt(w) for the pure Hubbard model 
for increasing interaction strengths: U/t = 2 (a), 3 (b), 4 (c), 6 (d), 8 (e), 12 (f), 20 (g). 

FIG. 3. For U/t = 12 and A/t = 3, (A) Disorder-averaged density of UHF single-particle states 
D(E) vs E/t (solid line), where the unperturbed (U = 0) bandwidth is 12t; IPR L'(E) for iV = 512 
site systems (dashed line); (B) Corresponding site-resolved local moment distribution |/x(e)| vs site 
energy e/t (solid line, left-hand scale), H(e;Ep) as defined in text (dashed line, right-hand scale). 

FIG. 4. Full transverse excitation spectrum for the Anderson-Hubbard model at RPA level 
(solid line) for N = 64 averaged over many disorder realizations, compared to the corresponding 
UHF Stoner spectrum (dashed) for fixed disorder A/t = 3 and U/t = 12 (a), 9 (b), 6 (c), 3 (d) 

FIG. 5. As Fig. | but with fixed interaction strength U/t = 12 and A/t = 1 (a), 2 (b), 3 (c), 
4.5 (d), 6 (e) and 7.5 (f). Note the persistence of a spin wave-like band for all disorder strengths. 

FIG. 6. Uniform RPA static transverse susceptibility Xu for fixed U/t = 12 as a function of 
disorder A/t. 

FIG. 7. Scatter diagram showing IPR L for low-frequency RPA excitations vs u/t for fixed 
interaction strength U/t = 12.0 and varying disorder, A/t = 1.5 (A), 3.0 (B), 4.5 (C), 6.0 (D). 

FIG. 8. Probability distribution of IPRs P(L) for U/t = 12 and A/t = 1.5 (solid), 3.0 (dashed), 
4.5 (dot-dashed) and 6.0 (dotted). 

FIG. 9. Probability distribution i7 RPA (e;L) (Eq. flO))) of the lowest N RPA excitations with 
IPR L over sites of energy e, averaged over excitations with L > 0.3 (solid line) and L < 0.3 (dotted 
line) for U/t = 12 and A/t = 3. 
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TABLES 

TABLE I. Spin wave character of 1ow-cj RPA excitations Q for the non-disordered and disor- 
dered Hubbard models. For A/t = 0, Q = r2(y max ); £1 is the average of Q,(u>) over the lowest N 



excitations. 

A/t = A/t = 3.0 U/t = 12 

u/t n u/t n A/t n 

1.0 0.02 3.0 0.05 0.0 1.00 

2.0 0.39 4.5 0.20 1.5 1.00 

3.0 0.71 6.0 0.33 3.0 0.98 

4.0 0.86 9.0 0.75 4.5 0.89 

> 5.0 1.00 12.0 0.98 6.0 0.88 
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